Quarkonium production in coherent hadron-hadron interactions at the LHC 
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The photoproduction of quarkonium in coherent hadron-hadron (pp/pA/AA) interactions for LHC 
energies is an important tool to investigate the QCD dynamics at high energies. In this paper we 
estimate the integrated cross section and rapidity distribution for J/^ and T production using the 
Color Glass Condensate (CGC) formalism. We predict large rates, implying that the experimental 
identification could be feasible at the LIfC. 
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I. INTRODUCTION 



In next future the Large Hadron CoUider (LHC) at 
CERN will start its experimental physics program. Cur- 
rently, there is a great expectation that LHC shall dis- 
cover the Higgs boson and whatever new physics beyond 
the Standard model that may accompany it, such as su- 
persymmetry or extra dimensions [l!]. However, several 
questions remain open in the Standard Model, which will 
be probed in a new kinematical regime at the LHC and 
determine the background for new physics. In partic- 
ular, the description of the high energy regime of the 
Quantum Chromodynamics (QCD) still is a subject of 
intense debate (For recent reviews see e.g. Ref. 01)- 
Theoretically, at high energies (small Bjorken-x) one ex- 
pects the transition of the regime described by the linear 
dynamics, where only the parton emissions are consid- 
ered, to a new regime where the physical process of re- 
combination of partons becomes important in the parton 
cascade and the evolution is given by a nonlinear evolu- 
tion equation. This regime is characterized by the limi- 
tation on the maximum phase-space parton density that 
can be reached in the hadron wavefunction (parton sat- 
uration), with the transition being specified by a typical 
scale, which is energy dependent and is called saturation 
scale Qsat ■ Experimentally, possible signals of parton 
saturation have already been observed both in ep deep 
inelastic scattering at HERA and in deuteron-gold col- 
lisions at RHIC (See, e.g. Ref. 0, Q). Although the 
geometrical scaling and the diffractive events observed 
at HERA, as well as the high-pT suppression observed at 
RHIC, have a natural interpretation in terms of the sat- 
uration physics, none of these phenomena can be taken 
as a conclusive evidence for a new regime of the QCD 
dynamics. This is due to the kinematical limitations of 
the experiments. Consequently, the observation of this 
new regime still needs confirmation and so there has bee 
an active search for new experimental signatures. 

In the last years we have proposed the analysis of co- 
herent hadron-hadron collisions as an alternative way to 
study the QCD dynamics at high energies [1, 0, [j, d, 
B [i3|- The basic idea in coherent hadron collisions is 



that the total cross section for a given process can be 
factorized in terms of the equivalent flux of photons of 
the hadron projectile and the photon-photon or photon- 
target production cross section. The main advantage of 
using colliding hadrons and nuclear beams for studying 
photon induced interactions is the high equivalent photon 
energies and luminosities that can be achieved at exist- 
ing and future accelerators (For a review see Ref. [ll|). 
Consequently, studies of -fp interactions at LHC could 
provide valuable information on the QCD dynamics at 
high energies. The photon-hadron interactions can be di- 
vided into exclusive and inclusive reactions. In the first 
case, a certain particle is produced while the target re- 
mains in the ground state (or is only internally excited). 
On the other hand, in inclusive interactions the parti- 
cle produced is accompanied by one or more particles 
from the breakup of the target. The typical examples 
of these processes are the exclusive vector meson produc- 
tion, described by the process jh Vh {V = p, J/^*, T), 
and the inclusive heavy quark production [jh XY 
{X — cc,bb)], respectively. Recently, we have discussed 
both processes considering pp 0, [1] , pA [§| and AA 0, 0] 
collisions as an alternative to constrain the QCD dy- 
namics at high energies (For recent reviews see Refs. 
[H, [ill)- Our results demonstrate that their detection 
is feasible at the LHC. However, some of these results 
were obtained considering the saturation model proposed 
by Golec-Biernat and Wiisthoff (GBW) several years ago 
[1J|. That model has been improved recently consider- 
ing the current state-of-the-art of saturation phy sics: the 
Color Glass Condensate (CGC) formalism [ll [H, [l^. 
This fact motivates a revision of some our previous es- 
timates. In particular, in this paper we revise our pre- 
dictions for J/^t production in coherent pp/pA/AA col- 
lisions. Another goal is to present, for the first time, 
predictions for the T production in coherent collisions 
considering the CGC formalism. 

This paper is organized as follows. In next section 
(Sectioninj we present a brief review of coherent hadron- 
hadron interactions, introducing the main formulae. In 
Section Hill we discuss the QCD dynamics and the satu- 
ration model used in the calculations. In Section Hvl we 
compare our predictions for the "fp cross section with the 
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HERA data and present the predictions for quarkonium 
production in hadron-hadron interactions. Moreover, we 
compare the current results to related approaches avail- 
able in the literature. Finally, in Section |V| we summarize 
our main results and conclusions. 
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where ^ is the differential cross section for the pro- 
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II. COHERENT HADRON-HADRON 
INTERACTIONS 

Lets consider the hadron-hadron interaction at large 
impact parameter {h > Rh^ + Rh^) and at ultra relativis- 
tic energies. In this regime we expect the electromag- 
netic interaction to be dominant. In heavy ion colliders, 
the heavy nuclei give rise to strong electromagnetic fields 
due to the coherent action of all protons in the nucleus, 
which can interact with each other. In a similar way, 
it also occurs when considering ultra relativistic protons 
in pp{p) colliders. The photon stemming from the elec- 
tromagnetic field of one of the two colliding hadrons can 
interact with one photon of the other hadron (two-photon 
process) or can interact directly with the other hadron 
(photon-hadron process). One has that the total cross 
section for a given process can be factorized in terms 
of the equivalent flux of photons of the hadron projectile 
and the photon-photon or photon-target production cross 
section In general, the cross sections for 7/1 interac- 



tions are two (or three!) order of magnitude larger than 
for 77 interactions (See e.g. Ref. Q). In what follows our 
main focus shall be in photon - hadron processes. Con- 
sidering the requirement that photoproduction is not ac- 
companied by hadronic interaction (ultra-peripheral col- 
lision) an analytic approximation for the equivalent pho- 
ton flux of a nuclei can be calculated, which is given by 
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where lu is the photon energy, 7l is the Lorentz boost 
of a single beam and ij = ujb/jL', Ko{r]) and Ki{ri) are 
the modified Bessel functions. Moreover, fj — ui {Rhi + 
RhJ/lL and Uifj) = Kl (fj) - Kl (77). The Eq. ^ will 
be used in our calculations of quarkonium production 
in pA and AA collisions. On the other hand, for proton- 
proton interactions, we assume that the photon spectrum 
is given by [TH . 
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with the notation = 1 + [ 

Q?ni„ = ^V[7i(i-2c^/x/^) 

The cross section for the photoproduction of a final 
state X in a coherent hadron-hadron collision is given 



(0.71 GeV^)/g 



(2) 



and 



'NN 



cess (7/1 — + XK)^ ujmi 
and V'S'nn is the c.m.s energy of the hadron-hadron sys- 
tem. Some comments are in order here. Firstly, the 
coherence condition limits the photon virtuality to very 
low values, which implies that for most purposes, they 
can be considered as real. Moreover, if we consider 
pp/pPb/PhPb collisions at LHC, the Lorentz factor is 
7l = 7455/4690/2930, giving the maximum c.m.s. 77V 
energy W^p w 8390/1500/950 GeV. Therefore, while 
studies of photoproduction at HERA are limited to 
photon-proton center of mass energies of about 200 GeV, 
photon-hadron interactions at LHC can reach one order 
of magnitude higher on energy. Consequently, studies 
of coherent interactions at the LHC could provide valu- 
able information on the QCD dynamics at high energies. 
Secondly, in pA interactions one have that due to the 
asymmetry in the collision, with the ion being likely the 
photon emitter, the photon direction is known, which 
will implicate an asymmetry in the rapidity distribution 
(see below). Finally, in this work we consider that the 
produced state X represents a quarkonium (J/^* or T). 
Since photon emission is coherent over the entire nucleus 
and the photon is colorless we expect that the events to 
be characterized by K = ft-i /12 and two rapidity gaps. 



III. QCD DYNAMICS AT HIGH ENERGIES 

Let us consider photon-hadron scattering in the dipole 
frame, in which most of the energy is carried by the 
hadron, while the photon has just enough energy to dis- 
sociate into a quark-antiquark pair before the scattering. 
In this representation the probing projectile fluctuates 
into a quark-antiquark pair (a dipole) with transverse 
separation r long after the interaction, which then scat- 
ters off the hadron ^195. In the dipole picture the ampli- 
tude for vector meson production reads as (See e.g. Refs. 
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A,(x,r,AX* ,(4) 



where ^E*^ ^(z, r) and 'I'J^^(z, r) are the light-cone wave- 
function of the photon and of the vector meson, respec- 
tively. The quark and antiquark helicities are labeled by 
n and n, variable r defines the relative transverse sepa- 
ration of the pair (dipole), z (1 — z) is the longitudinal 
momentum fractions of the quark (antiquark), A denotes 
the transverse momentum lost by the outgoing proton 
(t = — A^) and x is the Bjorken variable. Moreover, Aqq 
is the elementary amplitude for the scattering of a dipole 



of size r on the target. It is directly related to the S- 
matrix element S{x, r, b) and consequently to the QCD 
dynamics (see below). One has that [2l| 



Agq{x,r,A) = 



d'be- 



,6. A d<7q 



(5) 



where one have introduced the differential dipole-target 
cross section. Consequently, one can express the ampli- 
tude for the photoproduction of a vector meson in the 
final state as follows 

da, 



d^b 



(6) 



where T denotes the transverse polarization and one 
takes into account non-forward corrections to the wave 
functions [22,]. Finally, the differential cross section for 
vector meson photoproduction is given by 



da . 



1 
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^.(^^h^Vh) = —\A{x,A)\'il + p')^ 

where (3 is the ratio of real to imaginary parts of the scat- 
tering amplitude. For the case of heavy mesons, skewness 
corrections are quite important and the y a re also taken 
into account. (For details, see Refs. [20|.12]|'). 

The photon wavefunctions appearing in Eq. ([6]) are 
well known in literature [2l| . For the meson wavefunc- 
tion, we have considered the Gauss-LC model [21] which 
is a simplification of the DGKP wavefunctions. The mo- 
tivation for this choice is its simplicity and the fact that 
the results are not sensitive to a different model. In pho- 
toproduction, this leads only to an uncertainty of a few 
percents in overall normalization. We consider the quark 
masses rric = 1.4 GeV and mf, = 4.2 GeV. The param- 
eters for the meson wavefunction can be found in Ref. 
[2l| for the J/^* case. Accordingly, we have computed 
the parameters for the T case. 

The S'-matrix element S{x, r, b) contains all informa- 
tion about the target and the strong interaction physics. 
In the Color Glass Condensate (CGC) formalism [l5l.[T6l. 
[l7j . it encodes all the information about the non-linear 
and quantum effects in the hadron wave function. It can 
be obtained by solving an appropriate evolution equation 
in the rapidity y = ln(l/a;) and its main properties are: 
(a) for the interaction of a small dipole (r <C l/Qsat), 
S(r) « 1, which characterizes that this system is weakly 
interacting; (b) for a large dipole (r » 1/Qsat), the sys- 
tem is strongly absorbed which implies S{r) <C 1. This 
property is associate to the large density of saturated 
gluons in the hadron wave function. In our analysis we 
will consider the saturation model proposed in Ref. (2]| 
which generalizes the lancu-Itakura-Munier (IIM) model 
[23j . introducing the impact parameter dependence. In 
this model the differential dipole-proton cross section is 
parameterized by = 2Np{x,r,b), where 

. I ^ 

Np{x,r,b)- ) 



for f < 2 . 
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FIG. 1: (Color online) Energy dependence of the 7p cross 
section for quarkonium production. Data from HI and ZEUS 
collaborations JEEi. [j^ /. 



= 'rQsat{x,b) and the saturation scale is given by 
t{x,b) = (xo/a;)^/2[gxp(__^)]5^^ The expres- 



2-BcGC 

sion for f > 2 (saturation region) has the correct func- 
tional form, as obtained from the theory of the Color 
Glass Condensate (CGC) [l5|. Moreover, ^cs{x, r) — 
7sat + '"^y^^ is the effective anomalous dimension, which 

/bat ' KAy ' 

determines the behavior of the dipole cross section in the 
color transparency regime (f < 2) and introduces evolu- 
tion effects not included in the GBW model. The origi- 
nal parameters of the IIM model were reanalyzed in Ref. 
[2l| and are now given by: A = 0.159, — 5.95 x 10"'', 
No = 0.417, m„/d/^ = 0.14 GeV and = 1.4 GeV. 
Moreover, we assume in our calculations Bcgc — 5.5 
GeV-2. Hereafter, we label this model by b-CGC. The 
corresponding dipole-proton cross section is given by 
C!dip{x,r) = 2/ d%Np{x,r,b). 

For photon-nucleus interactions it is possible to extend 
the b-CGC model described above using the Glauber- 
Gribov formalism [24]. In this case, the nuclear scattering 
amplitude Na{x, r, b) will be parameterized as follows 



NA{x,r,b) = I 1 - cxp 



1 — exp [—a In^ [b r)] , for r > 2 , 



]^ATA{b) adip{x,r) 



where TA{b) is the nuclear profile function, which will be 
obtained from a 3-parameter Fermi distribution for the 
nuclear density. 



IV. RESULTS 

In Fig. [1] we present the predictions of the GBW and 
b-CGC models for the quarkonium photoproduction at 
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FIG. 2; (Color online) Rapidity distribution for J/'^ photo- 
production on coherent reactions at the LHC energy (see text). 



FIG. 3: (Color online) Rapidity distribution for T photopro- 
duction on coherent reactions at the LHC energy (see text). 



HERA and compare it with the current experimental 
data [2^, [2^. In particular, it is the first time that the 
b-CGC model is explicitly compared to the HERA data 
on vector meson production (J/vP and T). We have that 
in the HERA kinematical range both predictions are sim- 
ilar and describe the experimental data. However, these 
are distinct at larger energies, mainly in the J/^f case. 
This result motivates the reanalyzes of our previous pre- 
dictions for JJJ^ production in coherent hadron-hadron 
interactions [3, [9|. Moreover, the description of the T 
production implies that the saturation models are suit- 
able to be used in the calculations of this final state in 
coherent pp/pA/AA collisions. 

Lets calculate the rapidity distribution and total cross 
sections for quarkonium production in coherent hadron- 
hadron collisions. The distribution on rapidity y of the 
produced final state can be directly computed from Eq. 
Q, by using its relation with the photon energy w, i.e. 
y oc \ii{2Lu/mx). Explicitly, the rapidity distribution is 
written down as, 

dcr [hi + h2 ^ hi (g) X (g) h2] dN^{uj) 

1 — ^ J O'jh-^Xh [^) (oj 

ay auj 

where (g) represents the presence of a rapidity gap. Conse- 
quently, given the photon flux, the rapidity distribution 
is thus a direct measure of the photoproduction cross 
section for a given energy. In Fig. [2] and [3] we present 
respectively our predictions for J/^P and T production 
in coherent pp/pA/AA collisions, considering A — Pb 
and y^SNN — 14/8.8/5.5 TeV, respectively. In the pp 
case, the production at mid-rapidity at the LHC probes 
x-values of order (2 — 7) x 10^^. In the AA case, one gets 
(6-20) X 10-4. 

For pA collisions we have an asymmetric rapidity dis- 



tribution since the electromagnetic field surrounding the 
ion is very larger than the proton one, due to the co- 
herent action of all protons in the nucleus. As a conse- 
quence, the photon direction is known. In contrast, for 
pp and AA collisions we have symmetric rapidity distri- 
butions. In Table [J one presents the correspondent in- 
tegrated cross sections (event rates), using the expected 
luminosities. For comparison, we also present our predic- 
tions for quarkonium in the ultraperipheral collision of 
light nuclei. These results can be contrasted with those 
obtained for J/* production in Refs. 0, i, [13, [H, ill . 
Firstly, in comparison to our previous estimates for J/vt 
production in coherent pp and AA collisions , we have 
that the b-CGC predictions are almost identical to those 
obtained using the GBW model. It is directly associated 
to the fact that the main contribution for the total cross 
section comes from low energies, where the predictions of 
the GBW and b-CGC models for the quarkonium -fp pro- 
duction are very similar. In comparison to the estimates 
obtained in Ref. [2^ , our results give higher cross section 
by a factor of order 11% for Ca nucleus and are almost 
similar for Pb nucleus. On the other hand, in comparison 
with Ref. [2§| our results are almost 16% lower for Ca 
and 41% lower for Pb. As explained in detail in Ref. 
the difference between the predictions comes mainly from 
the distinct QCD approaches used and the different pho- 
ton flux considered in the calculations. Concerning the 
J/^f production in coherent pp collisions, we have that 
in comparison to that obtained in Ref. [27[ our predic- 
tion is almost 11% larger. Finally, for pA collisions, we 
have that our predictions are almost 12 % smaller than 
those obtained in Q , where we have used the IIM model 
[23| . which does not consider the impact parameter de- 
pendence of the scattering amplitude. 
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LHC 


J/* 


T 


pp (£ = 10^ mb-^s"^) 


136 nb (1360.0) 


0.8 nb (8.0) 


pPb (£ = 420 mb-^s"^) 


85.5 ^b (36.0) 


0.43 ^b (0.18) 


CaCa (£ = 43 mb-^s"^) 


440 fih (19.0) 


9.7 Atb (0.42) 


PbPb (£ = 0.42 mb-^s-i) 


43.5 mb (18.0) 


0.096 mb (0.04) 



TABLE L The integrated cross section (event rates/second) 
for the photoproduction of quarkonium in pp/pA/AA colli- 
sions at the LHC (see text). 



Lets consider now the T production. In this case, our 
estimate for pA are original in the Hterature using the 
CGC approach. In contrast, for pp and AA coUision it 
can be compared with those from Refs. [13, HO]. In 
comparison with Ref. [l^l we have that our results for 
pp (PbPb) are approximately a factor 4 (2) smaller. This 
large difference is directly associated to the very steep 
parameterization used in Ref. [2^ for the jp cross sec- 
tion [a-yp (X W]fp). On the other hand, our predictions 
for T production in PbPb collisions are approximately 
30 % smaller than those obtained in Ref. [30] consider- 
ing the impulse approximation. In comparison with the 
prediction obtained using the Glauber plus leading twist 
approximation, our result is approximately 20 % larger. 

Finally, two comments concerning the different models 
used to described the vector meson production and the 
hadroproduction of vector mesons are in order. Firstly, 
the vector meson production in general is described by 
quite distinct approaches for light and heavy mesons. 
In the light meson case, the old vector meson domi- 
nance (VMD) model is often adopted '28'|, whereas for 
the heavy mesons a leading logarithmic approximation 
of the coUinear approach is considered [29|. In this case, 
the hard QCD scale is given by the vector meson mass. 
On the other hand, the dipole picture and saturation 
physics approach used in this paper is the unique theo- 
retical formalism available that describes simultaneously 
light and heavy vector meson production [23|, with the 
transition between light and heavy mesons being dynam- 
ically introduced by parton saturation effects (via satu- 
ration scale) in the target. Secondly, the vector mesons 
can also be produced in inclusive and exclusive hadron 
- hadron interactions. In comparison with inclusive vec- 
tor meson production, which is characterized by the pro- 
cess hi + h2 V + X (For recent reviews see, e.g., 



|3ll.l33|). one have that the photoproduction cross section 
is smaller by approximately three order of magnitudes for 
proton-proton collisions. For nuclear collisions, this fac- 
tor is smaller due to presence of shadowing corrections. 
Although the photoproduction cross section would be a 
small factor of the hadronic cross section, the separation 
of this channel is feasible if we impose the presence of two 
rapidity gaps in the final state. This should eliminate al- 
most all of the hadroproduction events while retaining 
most of the photoproduction interactions. However, two 
rapidity gaps in the final state can also be generated in 
exclusive hadron - hadron interactions, characterized by 
the process hi+h2 hi(^V(^h2. Recently, the exclusive 
J/vP and T hadroproduction inpp/pp collisions were esti- 
mated in Ref. |33l| considering the pomeron-odderon fu- 
sion (See also [3J|). Although there is a large uncertainty 
in the predictions for pomeron-odderon fusion, it can be 
of order of our predictions for the photoproduction of 
vector mesons. As pointed in Ref. [33,] the separation of 
odderon and photon contributions ahould be feasible by 
the analysis of the outgoing momenta distribution. How- 
ever, this subject deserves more detailed studies, which 
we postpone for a future publication. 



V. SUMMARY 

The QCD dynamics at high energies is of utmost im- 
portance for building a realistic description of pp/pA/AA 
collisions at LHC. In this limit QCD evolution leads 
to a system with high giuon density. If such system 
there exists at high energies it can be proven in coherent 
pp/pA/AA collisions at LHC. In this paper we have an- 
alyzed the quarkonium production considering the most 
recent phenomenological saturation model and demon- 
strate that it describes the HERA data. We have inves- 
tigated the J/^* and T production in coherent hadron- 
hadron collisions and showed that the experimental iden- 
tification could be feasible. 
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